Theory of the tunnehng resonances of the bilayer electron systems in strong magnetic 
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We develop a theory for the anomalous interlayer conductance peaks observed in bilayer electron 
systems at v = 1. Our model shows the that the size of the peak at zero bias decreases rapidly with 
increasing in-plane magnetic field, but its location is unchanged. The I-V characteristic is linear 
at small voltages, in agreement with experimental observations. In addition we make quantitative 
predictions for how the inter-layer conductance peaks vary in position with in-plane magnetic field 
at high voltages. Finally, we predict novel bi-stable behavior at intermediate voltages. 
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PACS numbers: 73.43.Lp 

The discovery of dual peaks in the interlayer tunnel- 
ing conductance in double layer quantum Hall systems at 
total filling factor = 1 by Spielman et. al. M has stimu- 
lated a number of theoretical studies. ^, ^,]5|, |[ 0| Such 
interlayer tunneling measurements are a valuable tool to 
study the dynamical aspects of bilayer electron systems 
(BLES), since the incident quasiparticles interact with 
the tunneling barrier, impurities, interface roughness of 
the wells, as well as the two-dimensional electron sys- 
tem (2DES) in the wells. The last of these gives rise to 
an inelastic scattering mechanism where the tunneling 
quasiparticle excites the collective modes of the BLES in 
strong magnetic field. The dispersion relations of these 
modes, and hence the inelastic scattering rate are sensi- 
tive to the external in-plane field, producing a resonant 
peak in the conductance which varies in voltage with the 
applied in-plane magnetic field. The observation of this 
peak at non-zero voltage was been reported by Spielman 
et al. [0. In addition, tunneling quasiparticles can in- 
teract with topological defects in the order parameter 
such as merons (which carry the electrical charge) caus- 
ing phase decoherence and dissipation in the tunneling 
current with an inelastic scattering rate denoted by 
a±. The height and the width of the interlayer current 
peak is limited by such dissipative effects. 

The least understood aspect of the experiment is the 
yet unexplained peak in the conductance at zero voltage 
and its dependence on the in-plane magnetic field. The 
location of the zero bias voltage peak is insensitive to the 
in-plane magnetic field, but the height of this peak de- 
creases rapidly with in-plane magnetic field. Among the 
many theories of this experiment |[ ||, ||] there are 
different and even controversial interpretations for these 
observations. For example in Refs. [|[ ||, it has been 
argued that the conductance peak is the remnant of the 
long Josephson effect, however, in Ref. |6j the micro- 
scopic calculation indicates this system can be described 
by the excitonic superfluids. 

In this paper we propose an approach based on a 



damped Landau-Lifshitz equation for the pseudospin 
order parameter. Our model differs from other 
approaches Q in the types of dissipation included and 
how interlayer current is calculated. ^, ^, |6| These 
differences permit allow the model to capture previously 
unexplained features of the interlayer tunneling spec- 
trum in low and high bias voltages. Our model repro- 
duces the experimentally observed linear interlayer cur- 
rent {It ~ V) at low voltages (fig. 1), along with the value 
of the peak in the conductance at zero bias. Taking the 
intrinsic damping mechanism into account, we show the 
height of the conductance peak falls off as l/Q^, where 
Q = {B\\/ B±){d/£Q) is the in-plane magnetic field wave 
vector, £o is the magnetic length, and {B±) is the 
component of the the field parallel (perpendicular) to the 
quantum well. At higher bias voltages the system enters 
a non-linear regime with a conductance peak whose posi- 
tion changes with in-plane magnetic field. At still higher 
voltages, the current decays It Furthermore, we 

argue the low voltage state is different from the Joseph- 
son effect, and therefore argue against possibility of the 
Josephson effect at V — > 0. Moreover, we predict the ex- 
istence of a new bistable state between a "rotating" and 
"locked" states for the order parameter for small volt- 
ages in the presence of in-plane magnetic field, could be 
realized, depending upon the initial conditions. 

Our interpretation of the electrical current differs from 
other approaches. We model the steady state flow of 
quasiparticles, as an imperfect capacitor with a non- linear 
charging energy, hence the number of electrons and holes 
in different layers is fixed (but not equal). This system 
of a parallel resistor and capacitor (RC), connected to an 
external electrochemical potential gives the steady state 
dissipative interlayer current / = Sq{V)/Tz. Here 6q is 
the restored charge in the capacitor (which is a function 
of external potential V), and = RC is the relaxation 
time of the circuit. This analogy leads us to introduce a 
damping coefficient az{= l/r^) in the theory of the inter- 
layer tunneling effect. This coefficient, absent in previous 
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FIG. 1: The height of the zero voltage conductance peak 
(solid line) and the conductance peak at non-zero voltage 
(dashed line) versus the in-plane magnetic field, obtained by 
Eqs. (^ and (|^). Stars and circles are the experimental data 
for total electron density Nt = 6.0 x lO'^'^cm"^. Inset: the 
magnified crossing region of the main panel. The height of 
the zero- voltage peak falls off' as the non-zero voltage 

peak varies as 1/Q. 



models, is crucial to producing a zero bias peak. 

The low energy physics of the bilayer quantum Hall 
(pseudo) ferromagnets, which is based upon the micro- 
scopic Hartree-Fock model at odd integer filling fractions 
is described by the following effective Hamiltonian: 
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dx J \ dx 
^mxCos{Qx) + mySm^Qx)"^ , 




(1) 



where 'm{x,t) is the order parameter unit vector {ruz is 
the particle density difference between two layers, and 
rUx, and rUy are its canonical conjugate variables), pE is 
the in-plane (pseudo)spin stiffness, /3 gives a hard axis 
anisotropy due to the capacitance energy cost, Ag^s is 
the tunneling amplitude, V is the external interlayer bias 
voltage, and Q is defined above. 

Without the external leads, = is the lowest en- 
ergy state of an isolated BLES. Connecting the layers to 
the external leads bring this system out of equilibrium. 
Similar to an imperfect capacitor, the quasi-particles can 
flow between the layers, via a leakage current. Before 
reaching steady state, the displacement current drriz/dt, 
(which passes through the capacitor even at i^sAS = 0), 
can be measured. At steady state, the charge density 
on the capacitor becomes fixed, even though there is 
still a leakage current. The capacitance charge is given 
by < mz{x,t) > (<> is the average over the tem- 
poral and the spatial fluctuations), which is fixed and 
hence < dniz/dt >= 0. However, there is still a steady 
state current due to the interlayer quasiparticle tunneling 
channel. (This result is distinct from other approaches in 



which < drUz/dt >= would imply there is no tunneling 
current.) 

The energy loss by dissipative quasiparticle tunneling 
can be given phenomenologically by including a damping 
coefficient az, coupled to m^. 



It = azC < mz{x, t) > 



(2) 



The parameter which controls the resistance of the 
system is equivalte to the RC relaxation rate, and scales 
like l/y/(3AsAS- To calculate niz we start from the 
(iampec? Landau-Lifshitz equations (see e.g. fll]] ): 



dm 



^^cff — —2 



dt 

[ dn a_ dn 

1 dm dx 



(3) 



|, is the effective magnetic 



field, i? is a Landau-Lifshitz damping term (and h — 1). 
Our phenomenological model for the damping in BLES 
are characterized by two coefficients a^, and a± that 
define the damping vector R: 

Rx = —az'mlm'l — a±my{niym'^ — mr^niy), 
Ry = -azmlmy + a±mx:{mym'^ - rnxniy), 

Rz = azmz{mj;ml + myUi^). (4) 



which is defined so that rh ■ 
of rh is conserved (m • dtrh 



R — and thus the length 
- 0). The vector is the 
equilibrium value of the order parameter, ft is important 
to note the solutions of these Landau-Lifshitz equations 
exhibit different behavior as the bias voltage is increased. 
As — > (but AsAS 7^ 0) the tunneling term is dom- 
inant and the order parameter stays (almost) along the 
x-direction. Without damping the order parameter can 
precess around this direction, tracing out a cone centered 
on the TTix axis. The effect of the damping is to equili- 
brate the order parameter along the x-direction in a finite 
time, hence we assume m° = (1,0,0) in Eq.(^. Increas- 
ing the (small) bias voltage V alters the equilibrium state. 
Without damping, the lowest energy state can be deter- 
mined by minimizing Hamiltonian (|l|) . The role of damp- 
ing is to relax the excited states to some steady state m°. 
As V increases, the direction of ivP rotates towards the 
y-axis and tilts up slightly from xy-plane. The non-zero 
value of rriy reflects a non-zero Josephson current, but 
this current will vanish if V = due to the damping. 
We refer to this family of solutions as "damping-locked 
states." Starting from 7fi° = (1,0,0) in Eq.(^, and 
cranking up the bias voltages, when eV ~ ^/8f3AsAS, the 
electrostatic energy becomes comparable to the tunneling 
energy, and the order parameter starts to precess around 
a direction given by H^g. The system can no longer fol- 
low the tunneling term, and the order parameter becomes 
"unlocked" due to the bias voltage. The amplitude of 
these "unlocked" oscillatory solutions decreases with in- 
creasing V, so that at very large bias voltages the order 
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FIG. 2: Locations of non-zero voltage conductance peaks 
versus in-plane magnetic field. The theoretical curve (solid 
line) is derived from expression (^), while the filled circles 
are experimental data [7] for the same total electron density 
Nt = 6.0 X W^°cm-^. 

parameter aligns with the z-axis as 1^ — s- oo. In the rest of 
the paper we detail the solutions of the damped Landau- 



Lifshitz equations, and evaluate the interlayer current. 
Starting from the uniform solution mP, it is straightfor- 
ward to linearize the Landau-Lifshitz equations, to find 
their inhomogeneous solutions, using a starting point in 
our perturbative expansion that is different in low and 
high V limits (due to different nature of solutions). 
Large voltages: At high voltages {eV :s> \/8i3Asas) 
the pseudospin rotates around z-axis with a frequency 
uj = eV/h. It is then more convenient to work in a ro- 
tating frame which can be introduced by the transforma- 
tion n_L = m±exp{iu;t), Uz = rriz, and choosing n'^"^ = 
(1,0,0) which is equivalent to m^") = (coswi, sincji, 0) 
in the rest frame. Following this, the Hamiltonian (|l]) in 
the rotating frame can be transformed to ^: 

— As AS ^'n-xCos{ujt + Qx) + nySm{ujt + Qx)^. (5) 

Replacing fi" = (1,0,0) in Eq.(^), the Landau-Lifshitz 
equations can be derived: 



dfi ft 
a±nl - azul + = -ifSuyUz - 2p_E?^^-^^ - 2AsAsnzSm{ujt + Qx) 



a±nyn^ + = APn^Uz + 2pEnz-^^ + 2AsAsn^cos(wt -I- Qx) 



aznzTix + = 2pE < "x-Q^ ^ ^^~dx^ f ^ '^^SAs\^n^sm{ujt + Qx) - nyCos{u!t + Qx)j. (6) 



The last equation in Eqs. (|6|) can be interpreted as the 
continuity equation for the interlayer current. The exter- 
nal (but self-consistent) chemical potential contributes to 
the current via the first term in the left side of this equa- 
tion « 1). The first term in the right hand side gives 
the current density due to phase slips, J = pEdip{x)/dx, 
equivalent to a dissipationlcss supercurrent density of the 
excitonic condensation. In the presence of the small in- 
plane magnetic field (the commensurate state) J — peQ- 
Finally, the last term is analogous to the AC Josephson 
current. The perturbative solution around fp = (1,0,0) 
can be achieved by making the harmonic expansion: 
n = Asin[ijjt + Qx) + Bcos{ujt + Qx) + + ■ ■ ■. Substitut- 
ing this into equations (^, A and B can be determined 
after linearizing the Landau-Lifshitz equation, and one 
can derive the non-homogeneous leading terms in n(x, t). 
Plugging this into Eq. (||) (after replacing the coefficients 
A and B in h{x,t)), and making the space-time average, 
we finally end up with an expression for the steady state 



tunneling DC current 

r ^ SePAl^gujaz + a±) 

* {8f3pEQ^-Lo^ + aza^)^+u^az + ai_)^' ^' 

Solution (0) qualitatively well describes the peaks at 
UJ « ^/SPpeQ, corresponding to the resonance condition 
for the gapless acoustic mode. The height and the width 
of the interlayer current are controlled by the damping. 
In the absence of the damping /((oc A|^g) has a singular 
peak at ^/8/3peQ- Solution (^ is parametrically unstable 
for uj"^ > 8(3peQ^ ("tachionic" regime) as in the case of 
long Josephson junctions, (see e.g. Q), but it is stable 
for the large voltage limit lj^ ^ 8(3 peQ^ where the I-V 
characteristic follows the power law It ~ . Although 
it has been speculated by Fogler and Wilczek that the 
interlayer current peaks resemble the long AC Josephson 
effect [0 (where the location of the peaks are shifted by 
q;^, and ckz), here we argue the observation of these peaks 
is the manifestation of the spontaneous phase coherence, 
where the lowest energy state of the electrons is in a sym- 
metric linear combination of two layers which allow the 
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electrons tunnel through the energy barrier between two 
quantum wells without resistance (if a± = az = 0). It is 
also possible to search for the excitonic superfluid modes 
through out the interlayer tunneling measurement. The 
staggered supercurrent density of the excitonic pairs in 
the superfluid state is given by Jg = PeQs for low tun- 
neling energies. The velocity of the collective modes al- 
ters by Js : ujQ cjQ+Q^ . In the interlayer tunneling 
effect, the incident electrons are scattered by these col- 
lective modes. From the conservation of the energy and 
momentum, we have eV = Hujq+q^, i.e., the superfluid 
current shifts the location of the peaks. 

Close to the order-disorder transition point, the exper- 
iment suggests the possibility of the coexistence of the 
incompressible state with the compressible state ||l^. In 
this circumstances, the current from one layer to other 
layer can transfer through the phase coherent channel, 
as was described above, and also through the quasiparti- 
cle channels which are in the "uncorrelated" state. The 
contribution from the latter to the interlayer current in 
low voltages is dominant, which gives rise to a linear I- 
V characteristic, and therefore the possiblity of the DC 
Josephson effect is ruled out (see below). 

Small Voltages: For low bias voltages, we begin with 
the Landau-Lifshitz equations, Eqs.(||), and (^ (in the 
rest frame). Similar to Eqs.(^ we can derive a set of 
equations in this limit. The perturbative solution about 
m = (1, 0, 0) can be obtained by a calculation similar to 
that above. First we consider the simplest case of the zero 
in-plane magnetic field ((5 = 0). The uniform and static 
solution can be obtained easily. Given these solutions, 
one can find the interlayer conductance 



interlayer conductance can be obtained 



2e'^azAsAS 



4AsAs(2/3 + AsAs) + aza± 



(8) 



A similar technique can be used to derive the analyt- 
ical (non-uniform) solution in the presence of in-plane 
magnetic field if psQ^ ^ Asas- The effect of the tun- 
neling term (which is similar to the driving force) is to 
create only the harmonics with "wave-number" Q which 
itself will nonlinearly generate zero, 2Q and higher har- 
monics. Because of the damping, the amplitudes of all 
other harmonics (but the first harmonic) is expected 
to be zero in the zero temperature limit. It is there- 
fore natural to start with the following harmonic ex- 
pansion in the rest frame nix = 1 ^ where 
m+ = Ae'^^^ + Be~^^^ + + • • •, and m+ = mz+ iruy. 
Assuming that the leading perturbative terms should be 
the first harmonics of the driven wave-number we can de- 
termine A = -B = AsAsi4:(3 - ia±)/{8f3pEQ'^ +aza±). 
Substituting these into the Landau-Lifshitz equations, 
and linearizing them in terms of A, and i?, we can find 
the zeroth harmonic term tUq , and then the interlayer 
DC current It = a^e. < mz{x,t) >= azcm^o, and the 



Gt 



ie^PA^ 



SAS^ 



aza^{8ppEQ^ + azai_) + (32/32 + 2al)A 



-.(9) 



SAS 



The height of this peak falls off like (for high in- 

plane magnetic field), but the location of its center does 
not vary with in-plane magnetic field, consistent with |^ 
(see also Fig. 1). This perturbative solution is valid for 
small y's and large Q, and it coincides with the residual 
zero voltage peak in the presence of in-plane magnetic 
field. 

We note in passing that these solutions are valid in 
their respective limits, but that at intermediate voltages 
it may be possible to have more than one solution to a 
nonlinear differential equation. The basin of attraction of 
the solutions will depend upon damping and other details 
of the system. 

Numerical results: Our estimate shows the cross over 
between low and high bias voltages occurs at O.Olmy. 
To examine the accuracy of our model, the interlayer 
conductance peaks have been drawn, by using two ad- 
justable parameters. In Figs. 1 and 2 the fit to the ex- 
perimental data [Q is obtained by the following damping 
coefficients: a± — 0.250;^, a±az{a± + a^)'^ = 32f3'^A'gy^g 
(for AsAS = 90 pK, and ^/8/3Asas = 70mK, we find 
az = 75mK, and a± — 18mK). In Fig. 1, the height 
of the central residual (solid line) and the split off peaks 
(dashed line) vs. have been derived by Eq. (0), and 
(^). In Fig. 2 we present the locations of the split off 
peaks derived from Eq. (Q). 

We presented a physical picture based on a driven- 
damped easy-plane pseudospin ferromagnet model for 
the experimental observation of the interlayer conduc- 
tance peak in a bi-layer electron system at = 1 [Q. 
The first theoretical prediction for low bias voltage con- 
ductance peak vs. in-plane magnetic field has been made. 
It has been shown, at high voltages, due to the non-linear 
behavior of the capacitance energy, the inter-layer cur- 
rent shifts by in-plane magnetic field. 
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